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Real clocks are not perfect. This must have an effect in our predictions for the behaviour of a 
quantum system, an effect for which we present a unified description encompassing several previous 
proposals. We study the relevance of clock errors in the Zeno effect, and find that generically no 
Zeno effect can be present (in such a way that there is no contradiction with currently available 
experimental data). We further observe that, within the class of stochasticities in time addressed 
here, there is no modification in emission lineshapes. 
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I. INTRODUCTION 

' The problem of understanding time observables in quantum mechanics has a long and protracted history m . One of 
I the key observations in the process of a better formulation of time quantities was carried out by Misra and Sudarshan 
. 0, investigating how the measurement of Hfetimes could be affected by frequent probes into the evolution of the 
r"| ' system under study. It should be pointed out that Misra and Sudarshan placed their work in the context of time 
d ,: observables, making explicit connections to the problem of time of arrival, as discussed by Allcock 
' ' ' The result of Misra and Sudarshan that a continuously observed unstable particle would never decay was associated 
, by them to the name of Zeno of Elea, and it is under this title that the effect or paradox is currently known. In fact, 
CO . the effect had been pointed out before in different forms jH, 0] , as was indicated by Chiu, Sudarshan, and Misra . 
At any rate, the number of papers referring to the Zeno effect or paradox increased substantially after 1977, and even 
more so after the crucial experiment of Itano et al. (for a review of research in the topic up to 1997, with the 
corresponding bibliography, see Ref. |li|J ; a more up-to-date review of bibliography can be found in Ref. [ll|). 
\^ , Another intriguing aspect of time in quantum mechanics is related to decoherence and decoherence rates, both 
T-H in a general sense, and more specifically as a source of decoherence. Milburn |l2j | proposed a simple modification of 
^-H quantum dynamics in which the system does not evolve continuously under unitary evolution: it undergoes a sequence 
' of identical unitary transformations, which take place or not according to a Poisson distribution (i.e. the probability 
\ that there be n such transformations in a time interval t is given by a Poisson distribution). This proposal leads to 
' decoherence, while at the same time it conserves energy (a feature lacking in previous models of intrinsic decoherence, 
' such as the well known one of Ghirardi, Rimini and Weber jl^). 

A different aspect of stochasticity in time was put forward in Ref. |1J|: the fact that our clocks are not perfect 
-j-J ■ implies that incoherent superpositions of states at different instants of time are going to be necessary to account 
^ for the state observed at a given clock instant. The requirements for a clock to be considered good were examined 
^ ' and formalized, leading to a claim of uniqueness for the effective description of good clocks, in terms of a stationary, 
qh| Gaussian, and Markovian stochastic process. 

H^* ■ R. Bonifacio also proposed a generalization of Liouville's equation just by requiring that time be an stochastic 
variable and demanding that the a different form of stationarity hold 15j (see below, eq. JJl and the surrounding 
discussion, for the expression of this property) . He further posed a claim of uniqueness for the probability distribution 
' that encoded stochasticity in time, which he asserted was the F distribution. 

All these three proposals shared the result that, under adequate approximations, the systems behaving according 
to them would actually follow an evolution equation for the density matrix of the form 

p{t) = -^[H.p{t)]^^[H,[H,pm, (1) 

where we have set ?i to 1, as we do in the following, H is the Hamiltonian of the system under consideration, described 
with the density matrix p at clock time t, and k and § are constants with dimensions of time. In the formalism of 
Ref. jl3|, the quotient /■§ is associated to the intrinsic time step of the unitary evolutions; in Ref. 0| k stands 
for the strength of the correlation function of relative errors (the rates of increase or decrease of the clock error at 
different clock times), whereas •§ is the correlation time for those relative errors; to be complete we mention that the 
combination is equivalent to the "chronon" of Ref. [THj l. 

In fact, this kind of master equation was also known from the analysis of heat baths coupled to the s^tem by a term 
of the form HV , where H is the system's Hamiltonian, and F some bath operator (see for instance [igi section 2.3). 



The explicit connection between the proposal that quantum (and classical!) systems evolve according to non-ideal 
clocks, on the one hand, and the heat bath language, on the other, was shown in Ref. [T3 |. 

In a recent paper .17], Stephen Adler, using the language of Ito calculus, considered together both Zeno's effect and 
eq. QJ, with the result that the Zeno effect would be washed out by the new time scale jd. Given the results of 
Facchi and Pascazio (and others) 0, , it is rather surprising that the Zeno effect disappears no matter the value 
of the new time scale. Even more, if this equation is the result of a coupling between the system being considered 
and the heat bath, the total system (system plus bath) will of necessity present Zeno's effect. 

We should mention, for the sake of completeness, that the idea of errors in time measurements seems to be cropping 
up in several other contexts by various authors (for a couple of recent examples, see ^^ll^). The motivation behind 
most of those efforts seems to be either a desire to understand decoherence better, or the search for modifications of 
the ordinary quantum axioms regarding evolution. An example of this latter approach is to be found in the series 
of papers |2lL |23| ; notice that the specific form of non locality in time put forward in those papers leads to radical 
modification of the line shape, contrary to the results of Adler's and our own for non locality due to randomness in 
the measurement of time. 

The purpose of this paper is thus threefold: 1) to present a general formalism that accounts for various different 
proposals in a unified manner; 2) to examine whether Zeno's effect is indeed generically washed away by the mere 
fact that time should be considered as stochastic, or whether the result of Adler's does not indeed extend to more 
general situations; 3) to resolve the different claims of uniqueness which seem contradictory. 



A good starting point is given by the initial stages of Bonifacio's formulation which we present here with 
a notation and interpretation closer to that of Ref. [l^ . Let ps{s) represent the density matrix of system S at 
Schrodinger's ideal time s; that is to say, ps evolves according to von Neumann's equation 



where C is the Liouvillean (super)operator. 

This evolution equation is used to make predictions about the outcomes of experiments, that we can describe in 
the following manner: after a preparation stage, in which we make sure that we have set the system in a well defined 
initial state, we let the system evolve a time interval t that we measure with our clock. At that instant we measure 
the value of some observable of the system. Our predictions, in general, will not concern the actual value measured in 
a single instance of the experiment; they will rather provide us with probability distributions, which will be checked 
by many repetitions of the experimental procedure. 

However, any clock we might use will have intrinsic uncertainties, and the ideal time elapsed in each instance of the 
experiment will be different, even if we insist on always measuring the same time interval with our clock. We should 
then realize that the predictions we are required to provide must be predictions in clock time (<), not in ideal time 
(s), that is, we must have a way of computing probability distributions for all observables at time t. Through the 
standard arguments, this means that we can encode our predictions in a density matrix at clock time t, p(t). 

The actual observed/predicted density matrix at clock time t will be given by a superposition of the density matrices 
at different ideal times, under the assumption that, for different realizations of our experiment, the same reading t of 
our clock corresponds to different ideal time intervals, s, having elapsed since the preparation of the initial state. Let 
us describe this assumption by a probability density for those ideal times s at a fixed t, P(t, s): the actual density 
matrix, p(t), will be computed as 



Here one simple underlying assumption is that the preparation instant be labelled by s = and t = 0. Another 
assumption is that we know for certain that, at the instant we let the system go, the state is indeed always the same. 
Note that a simple way of modifying this assumption would be to allow s to be negative. When we introduce later 
licit J k), in eq. H10|l below, we shall actually relax this condition on s, as explained at that point. 

By the manner in which we have justified the introduction of P(i, s), we have already demanded that it be a 
probability distribution. However, an alternative way of arriving at the same concept would be to assume eq. (0) 
and that p{t) indeed be a density matrix (as is to be expected from the general arguments of Gleason's theorem — 
see the strengthened version provided by P. Busch ^^). It would follow that P(t, s) would have to be positive, and 



II. UNIFIED FORMULATION 



TPsi^) = -"^[H^Psis]] = -i£ps{s), 




(2) 



normalized to 1 when integrated over 5, 



/ dsP{t,s) 
Jo 



There is yet another important property that this probabihty density should fulfill: the way in which errors could 
accumulate should be independent of the instant of time at which we have started the clock. At least it should be so 
for good clocks; the way new errors are produced should be independent of the error up to that instant. Note however 
that the error at one clock time t + dt does depend on the error at time t. If this one is very big, it is very unlikely that 
at a later clock time the error could be zero. Therefore the stationarity requirement we are now discussing cannot be 
understood as stationarity for the errors in time measurement, but rather as stationarity in the buildup of errors. 

In order to obtain a mathematical statement of this stationarity requirement, consider the following setup: a system 
is prepared in an initial state ps(0), and is then evolved a clock time interval ti. The state at this clock time is now 



/>oo 

p{t,)= / dsiP(ii,si)e-'^^Ss(0). 
Jo 



By whichever means, the state is then frozen, and then we let it evolve a further time interval t2. This is another 
way of saying that we repeat the evolution, but now for a different time interval, and with an initial state that is now 
p{ti). The end result would be the density matrix 

p{ti+t2)^ ds2P{t2,S2)e-'''^p{h)^ ds2P{t2,S2)e~''-^ dsiP{h,si)e-'''^Ps{0). 
JO Jo Jo 

On the other hand, that freezing of the state need not be real, it is simply a tool of our imagination, and the density 
matrix at the end must be given by 



p{h+t2)^ / dsP{ti+t2,s)e~''^psiO). 
Jo 

This entails the following condition on the probability density P: 

P{ti+t2,s)= f ds' P{t2,S- s')Pih,s'), (3) 
JQ 

which is the mathematical expression of the stationarity requirement. 

In order to solve this functional equation, it is useful to use Fourier transforms. Define 



/•OO 

n{t,k) = / dsP{t,s)e' 
Jo 



Note for further use the following set of facts: since P is a probability distribution for the stochastic variable s, it 
follows that n(t, 0) = 1 (remember that s is stochastic with respect to the clock time which we measure); since P 
has support only on the positive half-line (as we are currently assuming) n(t, k) will have singularities only on the 
lower half of the complex k plane. If all moments in s of the probability distribution P existed, then Ii(t, k) would be 
analytic at fc = 0. 

By means of this definition, equation is transformed into the algebraic equation 

n(<i + t2,fc) = n(<i,fc)n(i2,fc), 

with solution 

n{t,k)^[g{k)f\ (4) 

where the dimensionful parameter A has been introduced to adimensionalize the exponent, and g(k) is an arbitrary 
function. Thus, A and g{k) fully determine the characteristic function Il{t,k), and consequently, the probability 
densities P{t,s). The physical significance of A and g{k) will become clear soon. For the time being, it is worth 
noticing that g{k) is the characteristic function of the probability density P{t, s) at time t — X, and that changes in 
the parameter A can be compensated by a change in the function g{k) as follows. Let g\{k) be an arbitrary function 
that determines n(i, fc) together with the parameter A, that is to say, under the condition n(A, fc) = g\{k). It is 
then straightforward to see that the same distribution n(t, fc) can be described in terms of a new function .9p(fc) and 
parameter fi, such that g\{k) — [^^(fc)]^^'^. 



The requirement that Tl{t, k) be the Fourier transform of a probabiUty density, with the definition above, imphes 
that g{0) = 1. Similarly, the condition on the location of the singularities also applies to g{k). Additionally, if we 
make the further hypothesis that all (s) moments of the density P(t, s) exist, this would entail that g{k) would have 
to be analytic at fc = 0. It is also clear that g{k) — > on the real line as |fc| oo, if P{t, s) is a continuous density. 

Assuming the existence of the first two moments (in s) of the distribution P{t,s), we can make the following 
interesting statements derived from the form of Il{t, k) expressed in eq. (@J. The expectation value of s is given by 



{s) = ds P{t, s)s = -idkUit, k) I 
Jo 



-*5'(0) 



where we have taken into account that ^(O) = 1. Observe that the expectation value is proportional to t, and that 
it is indeed exactly t if A = — ig'(O). It is immediate to appreciate that the ratio —ig'{0)/X measures the systematic 
drift of the expected values of time. The systematic drift can be eliminated by recalibration of the clock, identified 
by its characteristic function, and, indeed, the ratio —ig'{0)/X is invariant under the transformations X ^ fj, above, 
and pertains exclusively to the function n(t, k): if it is the case that A = —ig'^{0), it is then true that fi ~ — ig^(O). 
As to the variance, one can easily compute 

poo 1 

As2 = dsP{t,s){s^{s)f = j{-g"{0)+g'{Of)t, 

whence we see that As ^ \/i (with an adequate dimensionful proportionality constant). 

The simplest function that fits those criteria (to recap: ^(O) = 1, analyticity at A: = 0, all singularities in the lower 
half plane, and g{k) ^ as |fc| — s- cx) on the real line), and presents a singular point, is 

ffB(fc) = 



1 — ikr 

where t is a real number with dimensions of time. Notice that we are actually choosing a family of probability 
densities parametrized by A, i.e. by the instant of time at which the characteristic function Hit, k) equals gB^kl)- For 
each member of the family, as we shall see below, the meaning of the characteristic time r will be slightly different. 
The subscript B stands for Bonifacio: the probability density derived from this choice is 

1 p-s/r .„st/A-l 



TT{t/X) \t 

which is exactly the one put forward in Ref. The expectation value is (s) = Tt/X, and the dispersion As = T^JtjX. 
This means that among all the different members of this family of distributions, the one labelled by A = r is the only 
one for which there is no systematic drift in the expectation value of time; in which case r gives the rate of growth of 
the dispersion of the successive probability densities P{t^ s). 

Another simple alternative, with a very different analytic structure associated to periodicity, is given by 

5M(fc) = cxp(e''''^-l) , 

where, again, r is a real number with dimensions of time, and the subscript M now stands for Milburn. This function 
5m (fc) is an entire periodic function, bounded on the real line (but of course not everywhere, in keeping with Liouville's 
theorem). On performing the inverse Fourier transform, one is led to the discrete probability density 

oo / ^ \ " 



PM(t,s) = 2^^ h-) e-*/^5(s - nr) , 

which is precisely that proposed in Ref. (T2 | (eq. (2.7); in order to make the identification between this formulation 
and the original one it is convenient to use the construction of the next section). This discrete probability density 
has the interpretation provided in that paper: the probability that there be n identical unitary transformations of 
the density matrix, exp(— iT>C)/9(0) in a time interval t is given by the Poisson distribution. The expectation value (s) 
and the dispersion are given by the same expressions as those for (7b (fc)- 

The two g functions presented above do not of course exhaust all possible alternatives. Just as an example, consider 



where both t and a are positive time quantities (the subindex a stands for "alternative"). Note that if r and a are 
equal the resulting probability would be P^{2t^ s), which could also be phrased as -Pb(^, s) by halving the value of A. 
Without loss of generality, let t > a. The resulting probability density is 

with Iu{z) the modified Bessel function. This probability density serves as a counterexample to the uniqueness claim 
presented in Ref. . It should be pointed out that we make no claim whatsoever to the greater physical significance 
of this distribution, as compared to Pb , and they are to be evaluated according to their fitting whatever phenomena 
we would like to describe. The expectation value computes to be (s) = (cr + T)t/X, while the dispersion reads 
As = ^{t^ + a^)t/\. 



III. THE MASTER EQUATION 

Direct comparison of the second integral in equation (O with the definition of Il{t, k) tells us that we can write the 
averaging process leading to the observed density matrix at clock time t in the form 

p{t) - n(t, -C)ps{0) = e(*/^) '"[3(-^)Vs(0) . (5) 

By using the eigenoperators of the Liouvillean, that is to say, operators of the form |n)(m|, where \n) and \m) are 
eigenstates of the Hamiltonian, we can write the exact evolution of the components of the density matrix in that basis: 
let p{t) be given by ^^^^ ,^ Pnm\n){m\. The components of the observed density matrix at clock time t are related to 
those at time by 

P„m{t) =Tlit,Em - En)pnmiO) ■ (6) 

This formal exact solution is very useful when analyzing simple systems; however, it might be cumbersome in more 
involved situations, and some other simplifications and limits might come in handy. For this reason, let us rewrite eq. 
as a differential equation, 

m = \h,[gi-C)]p{t). (7) 
Since g(0) — 1 and g(k) must be analytic at fc = 0, eq. 10 admits an expansion whose first terms will be 

m - \ (-5(0)/: - i(-5"(0) + 5'(o)^)/:^ + . . .) Pit) , (8) 

or, taking into account that this expansion is actually the expansion of the generating function of the cumulants, 

m = (-^^c-\^c' + ..}jp{t). (9) 

In order for this expansion to be physically relevant it is clear that there must be either a renormalization of the 
energy, or else the first coefficient of the Taylor expansion of g{k) around A: = must be given by g'{0) = iX. In other 
words, the clock must be such that it adequately tracks ideal time, with no rescaling being necessary. Hence forward 
we shall assume that the statement ^'(0) = iX does indeed hold (which entails a relation among the parameter A and 
the characteristic times appearing in g{k)). It would be also pertinent to have some explanation for the validity of 
the expansion, which both for g^ik) and guik) reads, to second order in k, and making use of A = r (which is the 
specific form of the condition g'{Q) = iX), 

gB,M(fc) ^l + irk- T^k^ + O (r^fc^) ^ p(t) ^ _ l£2 ^ . . p(^t) . 

These examples lead us to understand the expansion as valid if the characteristic time of evolution (whose inverse 
gives us the characteristic expansion scale of £) is much larger than the time constants that appear in the definition 
of the g functions. In fact, those time constants, such as r [or a and r for ga{k)], characterize either the period (in 
cases such as guik)) or the closest singularities to the point fc = 0, so it is sensible to expect that this expansion will 



only be valid for characteristic evolution times larger than them, under the further demand that there is no resonant 
effect. 

In fact, Milburn provided an explicit example of the breaking down of the expansion due to a resonant effect, 
by examining the average value of an oscillator and noticing that some frequencies would lead to the freezing of the 
evolution of the oscillator, namely the harmonics of the evolution frequency 2it /t. Indeed, one should use in eq. ^ 



leading to freezing of the element Pnm{t) of the density matrix if Em — En = 27r//T, with I an entire number, i.e. in 
the presence of a resonance among the characteristic times of system and real clock. 



A completely different notation for the very same conceptual setup presented in section^]is that of path integrals. 
This will allow a more compact writing of some expressions, and will be of use in establishing a different set of 
approximations and connections with other models. The basic assumption is the following: a given clock time t can 
correspond to different ideal times s. These, therefore, give rise to a stochastic function of t, s{t), which we will 
find convenient to write as s{t) ^ t + A(t); that is to say, A{t) is the (stochastic) error associated with clock time 
t. The different probability densities P{t, s) lead to corresponding probability densities for the error variable. The 
stationarity requirement of eq. requires slightly more delicate handling. Let us go back to the use we have made 
of the probabilities P{t,s). We have set an initial time to, that we have so far denoted by 0, at which instant we 
are agreed that the elapsed ideal and clock time interval are both zero, and the density matrix is some initial /o(to) 
(which is for present purposes an unevolved matrix, and therefore there is no real distinction nor need to indicate the 
subscript S). Then some interval t has elapsed according to our clock, which reads the instant to +t; in ideal time an 
interval s has elapsed, associated to an instant in ideal time to + s. The stationarity condition means that we should 
not care whether Iq is 5 or 55, in order to determine the probability that an interval s in ideal time has elapsed when 
the clock interval is t. This condition can thus be translated into the following one on the error stochastic variable: 
that the relative error be a stationary stochastic process. In other words, that A(i) be written as /q rfCQ^(C): ^^'^ 
that this a{t) be stationary. Furthermore, since no negative time interval can ever elapse, we need to have a{t) > — 1 
always (this was also explained in 0| as due to "good" causality: if it were possible to have relative errors smaller 
than —1, we would see the clock going backwards) . We see that very often it is more convenient to deal with the 
relative error process, and we shall write our path integrals in terms of this stochastic function, with weight V[a]. 
In this manner, we recover the probability density P{t, s) as a sum over all possible sequences of relative errors (this 
is what J T^a will stand for), weighted with probability functional ^^[ck], restricted to those sequences for which it 
happens that t + A{t) equals s: 



Since a{t) > —1, the simplest possible choice for the functional weight V[a\, namely a Gaussian weight, is not 
feasible. Another way to see that this is indeed the case is to compute the Fourier transform of the probability density 




IV. ALTERNATIVE STOCHASTIC REALIZATIONS 





for a Gaussian weight, = exp [— ^ / d^d£^'a{^)c ^{^ — C')q^(C')] (foi' simplicity we have taken {a{t)) 

actually implies that the average value of s is t, thus choosing the correct scaling): 



0; this 



licit, k) 



I 



^ikt^-k^f{t) 



(10) 



where c(^) is the correlation function of relative errors. 




c ^(^) its functional inverse, and f{t) is explicitly defined in the last equality of eq. (|1(J|I . 



Condition can only hold for a Gaussian functional weight, as can be seen from comparison with eq. (|10() . if 
g{k) = exp(i/cA), which means that the clock is perfect and always gives perfect ideal time, i.e. P(t, s) = 6(t — s). 

On the other hand, a real clock could perfectly well go back in time, in which case the restriction above would no 
longer hold. The stationarity requirement ©, would need not hold either, since the imaginary process by which we 
arrived at it could now entail retrodictions. The generic form of the Fourier transform of the probability density for 
a Gaussian functional weight centered at (a{t)) = (as we must demand to avoid systematic slipping of the clock) is 
of the form 

as computed above, and the corresponding probability density 

where s is no longer restricted to being positive. The Gaussian process leads to an effective evolution equation 

PG{t) = (-iC- j{t)C')pG{t), 

which is precisely of the form of the master equation J^l, with a time dependent coefficient. The reason for the 
coefficient of the linear term in C to be —i is that the average value {a{t)) is zero. If it were not zero, then this 
coefficient would acquire a factor 1 + {a{t)). Concentrating now on the quadratic term in £, let us point out that 
the case of a constant coefficient corresponds to the Wiener process for the absolute errors. Notice that all the higher 
order cumulants vanish and the master equation is exactly given by the first two terms. 

This signals that whenever we can simulate the evolution of the quantum system by the master equation we 
can perform multi-time (and single time) computations using a Gaussian Markovian stationary functional weight 
7-" [a]: the Ornstein-Uhlenbeck process, in fact, in the stationary regime. Let us be more specific about this issue. If 
we indeed were to use the Ornstein-Uhlenbeck process for the relative errors, with correlation function {a{t)a{t')) — 
[n/'df exp(-|i - t'|/??), the function f{t) would be 



and the effective evolution equation would read 



Pou(i) = {-iC - ^ (l - e-*/'') Pou(i) , 



which indicates that there is a transient effect up to times of the order of the correlation time for the process, "d. 
Notice that at t = there are no non-Liouvillean terms, due to the specific character of the transient of the correlation 
function. As we shall see, this will be particularly relevant in Zeno's effect. After times of order i?, the evolution is 
dictated by the master equation Q). For the sake of later convenience, observe that the characteristic function of the 
master equation, under the identification r — 2k^/i5, is given by 

n,„.e.(t,fc) = e''^*-''-'^*/^ 
Quite obviously, the subscript m.e. stands for "master equation." 



V. THE ZENO EFFECT 



The generality of Zeno's effect has been shown in very many forms. An early clear description of the general 
character of the effect was provided by Chiu, Sudarshan, and Misra 8]: if the Hamiltonian is bounded from below, 
and the initial state is such that the expectation value of the Hamiltonian is finite, then the derivative of the survival 
probability with respect to time at the initial instant is zero, which entails that the decay of the survival probability 
is slower than any exponential. The survival probability at ideal time s is given by Tr{p(0) [exp(— is£)/9(0)]}. As, 
under these conditions, this survival probability has no linear term, close to s = 0, on performing sufficiently frequent 
measurements we find that the survival probability at any later time is equal to 1: the system is confined to its initial 
state. On the other hand, if the survival probability did indeed have a linear term in s, then the evolution under 



frequent measurements would be an exponential decay, with the decay constant given by the coefficient of the linear 
term. 

We shall now investigate whether the uncertainties in time do eliminate Zeno's effect from taking place. The 
quantity we want to investigate is the survival probability, that is, given an initial state p(0), we should compute 

p{t) = Tr [p{Q)p{t)] , 

where we measure the time elapsed with a real clock, and, therefore [see eq. (jSJ] 

p(t)=Tr{p(0) [n(t,-£V(0)]} . 

Note that the objects which include the Liouvillean are super-operators, not operators. This entails a slight compli- 
cation of notation, which we will fix by requiring that super-operators act on everything on their right. 

For measurements of time intervals that satisfy the stationarity constraint J^Jl, the linear term in clock time of the 
survival probability is given by 

jTv{piO)\n[gi-C)]piO)} , 
which, in terms of the energy basis, can be written as 

I |Pn™(0)|'ln|5(£:„-£;™)| . 

n>m 

Unless this quantity is zero, the quantum Zeno effect will not take place (notice by the way in the previous computation 
that, by construction, g{k)* — g{—k) for real k.). 

Let us assume the validity of the expansion in eq. lO, with (s) = t, and further that the initial state is a pure 
state, p(0) — \a){a\. In such a situation, the linear term will have a leading term of the form 

^^Tr[p{0)C'p{0)]=-^{AH)\ 

where {AHf = {a\H'^\ a) — (a\H\a)'^. That is to say, the small time survival probability will be 

A.52 

p{t) = 1 - — + ... , 
'''z 

with Tz = I /AH. No Zeno effect survives; frequent measurements of a system will not maintain it in the initial state, 
since As^ ^ t. There will be an exponential decay no matter how fast the measurements are! Admittedly with a 
very small decay constant, but exponential nonetheless. We can reobtain this result in a slightly different fashion, by 
noting that all terms of the Taylor expansion of the (ideal) survival probability in ideal time s contribute to a linear 
term in t when we perform the averaging with P(t, s): all cumulants are proportional to t, as we have seen, and this 
is itself a consequence of the stationarity property. 

For general clocks, that do not necessarily fulfill stationarity condition iPJ, the term linear in time will be of the 
form Tr [^(0)9(11(0, —C)p{0)], generically non vanishing. In the Gaussian case, which is not stationary in the sense of 
eq. Q, with initial pure state, leads to 

'^z 

Thus a very special situation is associated with the Ornstein-Uhlenbeck clock, for which /(O) ~ 0, and Zeno's effect 
survives. Note however the exceptional character of this case, due to the specific form of the transient. 

We have thus shown that clock errors generically wash out any possible Zeno effect. This being the case, how can 
we explain the experimental results of Itano et al. [p? Forgoing an analysis in terms of the full three-level system 
coupled to the electromagnetic field [2^ I25L l28| | . and neglecting optical pumping due to the measuring laser, we 

can concentrate on the following conceptual setup: consider a two level system undergoing a 7r-pulse Rabi oscillation. 
At regular intervals the system is queried as to the state it is in, whether the first or the second level. Formally, 
at regular intervals the coherences of the density matrix, p\2 and pi\ are set to zero. In the standard analysis, the 

evolution between measurements is ordinary unitary Schrodinger evolution with the Hamiltonian H ^ ^ 




Given an initial density matrix p(0) ^ bj' ^^'^ measurements at intervals Tr/{nQ,), the probability of finding 

the system in the second level f ^ j after time tt/H is 



If the evolution and measurement process were to take place according to a clock with characteristic function n(t, fc), 
the population of the second level would read 

1 ^^ .\ /TT 



P2 

where 

c{t)^^[n{t,n) + n{t,-n)] . 

It is worth mentioning that in the case of a perfect clock, for which npoifcct(i, k) — exp(ifci) both expressions coincide, 
as they should. 

In order to illustrate this result numerically, let us consider the characteristic function of the master equation, which 
leads to 

For the system examined by Itano et al., = 7r/256 ms^^, and r should be larger than 10^^ s for it to have any 
noticeable effect. As a matter of fact, r should be larger than 10^"* s for the effect to compete with optical pumping 
due to the measuring laser. On the other hand, the coherence times mentioned in % (550 s) lead to t < IQ-* s; 
additionally, the precisions mentioned for time quantities in their experiment are of the order 10~^ s or better. It 
follows that the inhibition of transitions from level 1 to level 2 due to frequent measurement is still in place, even 
though no "perfect" Zeno effect could take place if the clocks themselves are not perfect. 



VI. EXAMPLES: OSCILLATING AND DECAYING SYSTEMS 

Consider an oscillating system in initial state \a) evolving under the Hamiltonian H, 

H 

It is easy to compute the exact survival probability under unitary evolution with the preceding Hamiltonian, and then 
perform the averaging over clock errors, leading to 

^(^^ = ?F + T + 'f^'Re ^')] + ^Re [n(t, 217')] } , 
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where n' = yW+K^. 

The small time expansion of this exact evolution in the case of the decohering master equation with t = 2K^/'i9, 
is quite simply 

Pm.o.W ~ 1 - + 0(<2) . 

It is immediate to see that there will be no full Zeno effect, since there is a linear term in the expansion. 
If we were to use guik), the full result would be 



A particularly interesting aspect of this expression is that whenever fl' — 2m: /t the Zeno effect reappears: the motion 
of the system is indeed frozen, as Milbmn pointed out for different reasons. 
Finally, for good Ornstein-Uhlenbeck clocks we find 



pou(t) = JL (^i^4 ^ 1^4 ^ 2if2^2g-K^a-(t/^-l+e-'/'') ^ 1^4g-4K=a-(t/^-l+e-/-) ^os(21]'t)^ . (11) 

As soon as t is substantially bigger than the correlation time of the OU clock the two survival probabilities computed 
with the master equation and with the OU probability weight coincide. However, for small times their behaviour is 
radically different: there is no linear term in t in the expansion of po\j{t), in keeping with the general formal result 
presented above. 

Let us now perform an analogous computation for a model of a decaying system, in such a way that we can 
extrapolate without difhculty to more general decaying processes. On the by-side we will obtain an expression for the 
line shape, confirming generically the result of Adler's that it is not modified by phase decoherence 

Consider thus a system with a discrete orthogonal basis {|a)} U {|u;)}(^, where w takes values in some discrete set 
of frequencies, and such that the Hamiltonian can be written as H = Hq + V, with 

Ho = LOa\a){a\+J2^\^)H. (12) 

UJ 

and the only nonzero elements of V being (a|y|a;) and their complex conjugates. The initial state will be the pure 
state I a). This system is a simplified model of decay from this pure state to the rest of Hilbert space. 
By using the resolvent, G{z) = (z — H)~^, the exact Dyson-Schwinger's equation 

G{z) ^ Go{z) + Go{z)VG{z) 

for this system can be solved to 

CAE) = '""^l"' , ,,,, , (14) 

(E - UJ)[E -UJa- ^a(E)) 

where 

Ga{z) = {a\G{z)\a) , GJz) = {u:\G{z)\a) , 

and 

^'^(E)^^ (15) 

UJ 

is the exact self-energy for this model. In more general models, it is still the case that we can write the expectation 
value of the resolvent in state \a) in terms of the self-energy Ea(i?), as expressed in eq. (|13|) . The change will come 
about because of the modifications of the self-energy, which will no longer be determined by eq. (|15|l . The transition 
quantity Gi^ (E) will also have a different expression. 

The small time behaviour of the survival amplitude in ideal time s is determined by the large energy behaviour of 
Ga{E), and, consequently, of the large energy behaviour of the self-energy. In our specific example it is very easy to 
see that, for energies large in comparison to all to values. 



Sa(£;)«^El(«l^l^)l' = ^(H^' 



2ia) - (a|i/|a)2^ - ^ 



Under general assumptions, and for large energies, the self-energy function has exactly the same structure Y.a{E) ^ 
1/t'^E for generic models. Then, since the survival probability amplitude A{s) is determined by the resolvent and 
hence by the self-energy through 



A'^) = :^ I dEe~'^'Ga{E), 



where by B here we denote the adequate integration path, the survival probabihty at small ideal time s is 



Pidoal(s) = \A{s)\" 



4 + w2r| 



TZ 



(16) 



On performing the averages for Eq. (|16|l . we obtain 

Pit) - 



4 + 



Ren(t, .\/c.2+4/r|) 



The objection might be posed that small t does not imply that s is small (the comparison term in order to state 
the smallness or otherwise of these dimensionful quantities is always taken to be the characteristic evolution time of 
the system). However, we have seen above that clocks that satisfy the stationarity requirement of eq. display 
an average value of s that tracks t, and that the dispersion is given as As^ t. It then follows that indeed we can 
approximate small t behaviour by extracting the small s behaviour and then averaging. 

We see again that through this method we recover again that there is a linear term, and that the disappearance of 
Zeno's effect is completely generic. An exception can be found in the OU clock case, for which 



Pov{t) 



1-72. 



as we already know. 

We shall now study the large time behaviour of the system. Assume now that the states orthogonal to the initial 
one, |a), form a continuum, in the manner postulated by Weisskopf and Wigner |29l |. so that an imaginary part can 
arise for the poles of Ga{E). Then, for large times, the dominant behaviour, both for the survival amplitude and for 
the probability amplitude to find the system in a state different from the initial one, will be determined by the pole of 
Ga{E) closest to the real axis. Assume there is only one such relevant pole, of the form ujp — i'y/2. The contribution 
of this simple pole to the survival amplitude is 

Apis) = ^e-'^^'-^'/^ , 

where Zp is the relevant residue. Hence the contribution of the simple pole to the survival probability at ideal time 
s becomes 



Pp,idcal(s) = -2pe '^'^ . 

If this were the only contribution to the survival probability, it would be immediate to conclude that at clock time t 
the survival probability would read 

Pp(t) = ZpH(i,z7). 

For large times it is to be expected that the main contribution of real clocks will be given by the master equation 
approximation. It results that 



Ppit) « Zpc' 



'7(l-7T/2)t 



which means that the long times delay is slowed down as an effect of clock errors, the new decay constant being given 
by 7(1 — 7t/2). As the half-life and the error dispersion parameter approach, the long time decay will become slower! 

Under the same approximation (single simple pole or Weisskopf- Wigner's approximation), we have for the proba- 
bility amplitude of states orthogonal to the initial one, at very late times, and for the specific Hamiltonian (|12f) . 



BUs) = {L^e-''"\a) 



{uj\V\a) (e-^'^p'e 



LUp — LU — ij/2 



whence it follows that the probability of finding state \uj) at (sufficiently large) clock time t is 



Puit) = \BUs)\' 
= dis) 



MV\a)\' 



{ujp ~ ujX + 72/4 



(1 + H(t, ^7) - H(t, i^^ujp^ i-i/2) - n(<, ujp-uj + i-i/2)) 



MV\a)Y 



(wp- 6^)^+72/4 



which means that the hne shape, within the approximation carried out (Weisskopf- Winger approximation, see ref. 
p9j ) does not change at all because of the errors in the clocks. We have not actually fully proved this statement: 
we have not justified that the use of the long-time approximation in ideal time is enough before taking the average. 
However, for probability densities that lead to condition the average value of s is proportional to t, and its 
quadratic dispersion goes instead with ^/t. It follows that one can carry over the approximation of large ideal times 
to large Schrodinger time. 

VII. CONCLUSIONS 

We have shown using a general description of errors in clocks, 1) that decoherence, although a slow process if the 
clocks are good, starts taking place immediately for generic models of clocks; 2) that, as a consequence, Zeno's effect 
would never fully freeze the system in its initial state, for generic clocks. 3) We have found two kinds of exception to 
this result: whenever the decoherence itself freezes the state, as happens for some frequencies in the case put forward 
by Milburn, and in the Ornstcin-Uhlenbeck clock. 4) That, under the assumptions of real clocks put forward in 
and in this paper, there is no change in the line shape due to non locality in time. 
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